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Abstract 

Starting from determinants at finite temperature obeying an intermediate bound- 
ary condition between the periodic (bosonic) and antiperiodic (fermionic) cases, we 
find results which can be mapped onto the ones obtained from anyons for the second 
virial coefficient. Using this approach, we calculate the corresponding higher virial 
coefficients and compare them with the results known in the literature. 
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In two dimensional space one finds fractional statistics which depends on a continuous 
parameter defined on an interval where bosons correspond to one end and fermions to the 
other |]-@]. One of the most striking properties of fractional statistics is the prediction 
of fractional charge, i. e., noninteger multiples of the fundamental electron charge e, at 
low energies inside condensed matter differently from quarks which can only be probed 
at very high energies and are confined to hadrons. In particular, Laughlin explained the 
essentials of the fractional quantum Hall effect (FQHE) || |10[ assuming the existence of 
currents carrying fractional charges e/3 inside two dimensional samples, which has been 
confirmed in very recent experiments [TT|, [TJ| . 

Different formulations of fractional statistics were proposed, among which the most 
fundamental ones are the anyon formulation based on a quantum mechanical model con- 
sisting of charged flux tubes confined to move on a plane || and the axiomatic Haldane's 
fractional exclusion statistics which is based on very general assumptions on quantum 
statistical mechanics in arbitrary dimensions |13| . However, these two formulations are 
not equivalent in general, as it was soon noted by Haldane. One way to see this is to 
make explicitly calculations, e. g., for the virial coefficients on both formulations and then 
compare the results. In the case of the exclusion statistics it was shown that for a gas 
in a low-density expansion only the second virial coefficient depends on the correspond- 
ing statistical parameter while the others are exactly equal to the bosonic and fermionic 
values [|TJ], [OJ. On the other side, in the charged flux tube or anyon formulation, the 
second virial coefficient was found exactly || but the higher virial coefficients are not 
exactly tractable, apart form the fact that the third virial coefficient has an exact mirror 
symmetry in its dependence on the statistical parameter in respect to its semi-interval 
W, P7| . The calculation of higher virial coefficients for anyons have been attacked with 



the aid of perturbation theory and numerical methods showing that in general, there are 



corrections with respect to the corresponding bosonic or fermionic values [|I"8| -|27|1. 

Recently, we introduced a generalized partition function interpolating that ones for free 
relativistic charged bosons and fermions which we showed to lead, in the nonrelativistic 
limit, to results which are mapped onto the conventional ones for anyons, in particular 
for the second virial coefficient [p8|| . In this letter we are going to use this approach to 
calculate the corresponding higher virial coefficients. We find that they are all dependent 
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on the interpolating parameters, so they only coincide with the results from the exclusion 
statistics for the bosonic and fermionic cases. However, these coefficients are remarkably 
close to the ones obtained from anyon calculations. 

The route that we are going to take here follows Ref. p8] so that we use few but 



essential information which are necessary for a thermodynamical description of a system 
whose properties interpolate those corresponding to fermionic and bosonic gases. These 
informations are: first the (relativistic) particles should obey energy conservation and so 
despite they are described by some involved specific (say, 2+1 dimensional) equation they 
should also obey the Klein-Gordon equation. Hence, the corresponding partition func- 
tion can be written in terms of the determinant of this operator; second, this determinant 
should be calculated using some appropriate boundary condition; and finally the partition 
function is defined by a certain power of this determinant. We choose a boundary condi- 
tion that continously interpolates the bosonic and fermionic cases. Such a condition, in 
the Euclidean time compactified to the interval [0, /3], (3 being the inverse of temperature, 
can be written as 

#;^e^;f), (1) 

where 9 is a real continuous parameter. The spatial analog of this condition was discussed 
before in the study of the vacuum structure in gauge theories [^9) and to describe the 



spatial behavior of anyons in 2+1 dimensions, which acquire an arbitrary phase when 
turned around each other [p], H |7|. In a different context, a similar condition had been 
used to interpolate continuously between the Euclidean and Minkowski metrics to discuss 
a possible dynamical origin for Lorentzian spacetime |3(| [H]]. Note that when 9 = 
the above condition reduces to a periodic (bosonic) boundary condition and if 9 = 7r it 
becomes an antiperiodic (fermionic) one. 

So, we can write a partition function according to these prescriptions as (h = c = 1) 

Z = det a (-D 2 + M 2 ) e , (2) 

where — D 2 + M 2 is the Klein-Gordon operator and the covariant derivative is defined 
as D v = (do + ifi,di), with a shift in its time component due to a nonzero chemical 
potential /i introduced to describe a finite dinsity of (noninteracting) charged particles of 
mass M in (iV + 1) space-time dimensions. The subscript 9 means that the eigenf unctions 
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of the Klein-Gordon operator satisfy the boundary condition ([!]). The power o of the 
determinant enables us to correctly describe the partition functions for charged bosonic 
(9 — 0, a — —1) and fermionic {9 = ir, a = +1) gases. 

The determinant (|2|) can be exactly calculated, so that the corresponding free energy, 
Q(/3,fi) = — In Z, can be written as 



Sl(P, n) = 2Va — - cos(n9) (e n ^ + e"^) - Kn+i (npM) , (3) 

V 27r P/ n=l 2 

where K v (x) is the modified Bessel function of the second kind. If we particularize the 
parameters a and 9 to the bosonic and fermionic cases we shall find precisely the results 
known in the literature p2|-p5|1. 

Let us now restrict our analysis to the particular case of most interest which is the 
two dimensional (N = 2) case in the nonrelativistic limit since this is the context where 
anyons and their applications appear. The nonrelativistic limit in this situation is taken 
by imposing the low temperature limit (3M » 1, since we are restricting the allowed 
energies to be much less than the rest energy of the particles. Further, we consider only 
the particle content of the above free energy, corresponding to the fugacity z = exp 
So, disregarding the antiparticle content with z' = exp (— we find that 

V 00 

^(/M = --rEWK\ (4) 

P n=l 

where b n (/3) are the cluster coefficients: 

= " cos(n0) exp(-nPM) ; {(3M » 1). (5) 

Now, using the standard formulae which relate the virial coefficients a n with the cluster 
ones 



36|, M ( a i = 1 

He = 



a 2 = - TTX2 ; ( 6 ) 
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«3 = -2 7r ^- +4^; (7) 



(&i) 3 (6i 
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and substituting (EJ) into these relations, we find the virial coefficients corresponding to 
the generalized partition function (|2|): 

\2 

a 2 = — (1 -tan 2 0) ; (11) 

\4 

a 3 = (1 + 6 tan 2 + 9 tan 4 0) ; (12) 

36a 2 

\6 

a 4 = -—^ (tan 2 + 6 tan 4 + 5 tan 6 0) ; (13) 

A 8 

a 5 = j (1 - 60 tan 2 0-1170 tan 4 0- 2700 tan 6 0- 1575 tan 8 0) ; (14) 

3600a 4 v ; ' v 7 

A 16 1 5 145 35 21 

a 6 = r ( tan 2 H tan 4 8 H tan 6 H tan 8 H tan 10 0) , (15) 

6 a 5 v 288 24 144 24 32 1 ' v 1 



where A = ^2ir(3/M is the thermal wavelength in natural units. Higher order coefficients 
can be analogously calculated, but for the moment the above ones are sufficient. 

Before we compare these coefficients with those corresponding to the anyonic formu- 
lation, let us first note that a relation between the parameters 9 and a can be found, 
based on the fact that we do not expect to see any statistical behavior in the one particle 
partition function Z\ = Vb%. Once b\ is given by (ED with n — 1, we have to impose that 
cr cos 9 = constant in order to satisfy that condition. Then, to reproduce correctly the 
bosonic and fermionic results we choose 

a(9) = -cos' 1 9. (16) 



Using this relation we can rewrite the virial coefficients (|TT])-(15) for instance, as functions 
of a = a{6): 

a 2 = ^-(--a); (17) 
4 a 

a 3 = ^(4- 12 + 9(r2 ); ( 18 ) 

36 o~ z 

a4=^(-- + 9<7-5<7 3 ); (19) 
16 cr 
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^ A8 (-22b-i^ + i^ 2+ ir 4 > ; < 20 > 
*= A " < 21 > 

So, from eqs. (|lT|)-(|15|) or eqs. (p~7| ) - (|21| ) we clearly see that these coefficients are singular 
for 9 = (2k + l)7r/2 with k integer or equivalently when a = 0. This is not the case 
of the anyon higher virial coefficients which have been shown using perturbation and 
numerical methods to be finite to all orders except for the second where cusps are present 
at the bosonic points ]T7|-p7|. However, our coefficients and the usual anyon ones share 



other interesting properties. First of all, they naturally reproduce the usual bosonic and 
fermionic virial coefficients when we put 9 = and ir, respectively, so that 

a 2 F B = ±-\ 2 - a 3 = -U 4 ; «4 = 0; a 5 = -^A 8 ; a 6 = 0; . . . (22) 
4 36 3600 

and all coefficients (pTD-(pT5|) or equivalently (|I~7|) -([21~]) are periodic functions of 9, as ex- 
pected. Further, when 9 is away from the singularities, i. e., near the bosonic 9 = 2kn 
and fermionic 9 = (2k + 1)tt points these expressions give good approximations to the 
higher virial coefficients of anyons in a very simple way, as we are going to show bellow. 
Near these points we can use the approximation 

tan 2 9 ss sin 2 9 (23) 

so that we have the approximate coefficients 

A 2 

a 2 ~±— (l-sin 2 #) ; (24) 
A 4 

a 3 ~ — (l + 6sin 2 # + 9sin 4 fl) ; (25) 
36 

A 6 

a 4 ~ =F — (sin 2 9 + 6 sin 4 9 + 5 sin 6 9) ; (26) 
16 

A 8 

a 5 ~ (1 - 60 sin 2 0-1170 sin 4 9 - 2700 sin 6 9 - 1575 sin 8 9) ; (27) 

5 3600 v ) , K ) 

a 6 ~ T \ 10 (— sin 2 9 + — sin 4 9 + — sin 6 9 + — sin 8 9 + — sin 10 9) , (28) 
6 v 288 24 144 24 32 >> K J 

where the upper signs correspond to fermions and lower to bosons. Note that the odd 

coefficients have the same global sign for bosons and fermions according to Eqs. (|lT1)-(|T5|). 
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To stablish the relations between these coefficients and those that come from anyons, 
let us first write the exact result of Arovas, Schrieffer, Wilczek and Zee [||] for the anyon 
second virial coefficient: 



Ao 



X 2 



-2±2)|d -2a z 



(29) 



where upper signs correspond to fermion, lower to boson based anyons and a is the 
statistical parameter, which is assumed to be a periodic function of another parameter, 
say an angle 9'. 



Higher anyon virial coefficients have been investigated by several authors |17|-|27|1 with 
different approaches. In particular, Dasnieres de Veigy and Ouvry |^2, ^| using a highly 
involved perturbation theory were able to find corrections of order a 2 for the anyon 
coefficients up to the sixth one: 

„2~ 



A» = A 



1 a 
36 + 12 



(30) 
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(31) 
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QVQ 3-VQ QVZ y/3-l 



(32) 



A e , 



\ 10 a 2 



+ 



5 , y/2 + 1 
In 



65 , V3 + 1 
In 



In 



576 ' 864^ y/2 - 1 ' 288^ '"^-1 ' 27 V$ "'v^-l 

25 , 3 + 

In t= + 



9 , 4 + V10 
m 



A 10 a 2 I =F 



48^/6 3-^6 32^ 4 - V^O 



(33) 

576 1027 v ; 

following the same notation used in Eq. Using very precise Monte Carlo calculations 

these results were later corrected to order a 4 for the third virial coefficient by Maskevich, 



Myrheim and Olaussem |26| 

A, = A 4 



sin 4 9' 



1 sin 2 9' 
36 + 12vr 2 ~ (621 ± 5)vr 4 



(34) 
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where sin^' = ir\a\, and more recently by Kristofersen et. al. to the fourth anyon virial 



coefficient 27 



A 4 = X b 



^-f (4= ln (v^ + 2) + cos 9') + sin 4 0' (c 4 + rf 4 cos 5') 
16n z \y/3 J 



(35) 



where c 4 = -0.0053 ± 0.0003 and d A = -0.0048 ± 0.0009. 

Note that there is a small difference between the expressions of Aa, Eqs. ([Jl|) and (|35|), 
in respect to the coefficient of term sin 2 9' because of the substitution of ^1 by cos#'. This 
does not modify our analysis since we are considering 9' near 7r or zero. More generally, 
it seems to be a basic difference between these two approaches to higher anyon virial 
coefficients. In numerical simulations one usually works with single smooth functions 
which do not distinguish the statistics, while in the perturbation theory, as expressed in 
coefficents Aa and A e , Eqs. (|31|) and (|33|) , these coefficients are statistics dependent as it 
happens to the well known Arovas et. al. result, Eq. (|29|). 

So, near the fermionic point in our formulation (9 m n) and imposing the equality 
between (p4j) and ( p9|) with fermion signs, we obtain 



sin 2 9 ~ 2a 2 ~ — sin 2 & (36) 

IT 2 



which shows that our virial coefficients , Eqs. (p5|)-(^q), are very close to the ones coming 
from anyons, Eqs. fl30l)-(35). An analogous situation occurs if we start near the bosonic 



point (9 fa 0) and we identify 



sin 2 ^ ~4|a| -2a 2 ~ — | sin 6> / 1 - 4r sin 2 & (37) 



so that similar results are obtained but since the above relation include linear terms in 
| sin#'| one would also find odd powers of this term in the analysis which did not occur 
for the coefficients near the fermionic point. 

These results are in good qualitative agreement with the perturbative and numerical 
ones for anyons but quantitatively the numerical factors in our coefficients are greater 



than those appearing in Refs. [f22| , |23| , f2q , p7| . However, it is important to emphasize that 



the our coefficients have the expected properties like the exact mirror symmetry |T7j and 
also contain the particular fermionic and bosonic cases when we choose appropriately the 
parameter 9. 
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Another good test for our coefficients is to compare them with the ones coming from 
the average field formulation [j37H 39lj8ll where anyons are supposed to be a large collection 
of noninteracting bosons [y = 9' /n) or fermions [y = 6'/tt — 1) of density p in an average 
magnetic field, B = vp&o, where $o = h/e is the elementary quantum flux. In this case, 
the virial expansion reads 

^-, ± ^ + ^V^V--i^V + ---. (38) 



Our results (p4j)-(p8|) are remarkably close to these ones, suggesting that our results may 
be an intermediary formulation between the conventional ones for anyons and the average 
field approach. 

We can go even further in the calculation of higher virial coefficients within our ap- 
proach since we already found the cluster coefficients exactly. For simplicity, we just 
present the basic behavior of the n-th virial coefficient: 

b n ~ l 



n 7T 



^+e<iik-p > (39) 

i 3 

where c" and df are well known constants that can be determined from standard Eqs. 
(P)-(|T0|) and so on. Using eq. (|5|) one finds 



In 



^^tan 2i - 2 0, (40) 
i=i 

where the constants can be determined from Eqs. (|5|)-(]TD|)..., combined with (|^). 
In particular, (— l/2n) ra_1 fc™ are identified with the Bernoulli numbers. Note that these 
polynomials in tan 2 9 with even n > 2 vanish for 8 = and tt, reproducing correctly the 
usual bosonic and fermionic virial coefficients in two dimensions. 

It is the hope of the authors that the higher virial coefficients presented here, which 
interpolate the bosonic and fermionic cases, can perhaps shed some light on the thermo- 
dynamics of fractional statistics. It should be remarked that we discussed here the free 
case only. However, in more realistic models one should take into account the coupling 
between the particles and/or the interaction with external fields. Another point of interest 
is the possibility of connection of these results with a more fundamental theory starting 
from a Lagrangian density with bosonic or fermionic fields coupled to a Chern- Simons 
term at finite temperature to see if the quasiperiodic conditions would come directly from 
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the gauge symmetry of the problem. The discussion of these points are presently under 
investigation and will be reported elsewhere. Let us also mention that our approach is 
independent of the space-time dimension, as it happens in Haldane's exclusion statis- 
tics, although here we discussed the particular case of N = 2 in more detail to stress its 
relations with the anyon formulation on the basis of the virial expansion. Within this 
context, it might be interesting to discuss the relation of our approach with systems in 
other dimensions, e. g., the one-dimensional cases presenting exotic statistics. 

H.B.-F. thanks D. Bazeia for reading a preliminary version of the manuscript and for 
useful suggestions on it, and the hospitality of the Center for Theoretical Physics, MIT 
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by CNPq (Brazilian agency). 
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